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SUMMARY 

^ This  example,  which  is  one  of  the  simplest  to  exhibit  vortex  refraction  of 
sound,  is  solved  using  an  approximation  which  is  uniformly  valid  at  all  frequencies 
and  small  Mach  numbers.  The  region  of  solid  body  rotation  is  bounded  by  a vortex 
sheet  which  implies  that  the  results  obtained  have  little  physical  relevance. 
However  this  flow  field  gives  acoustic  equations  which  are  comparatively  easy  to 
solve  and  thus  an  ideal  test  case  for  the  various  approximate  techniques  used  to 
investigate  the  effect  of  more  realistic  vortices  on  sound. 

The  present  Report  gives  an  assessment  of  a method  which  sets  out  to  express 
the  acoustic  pressure  as  an  asymptotic  power  series  in  the  Mach  number.  It  is 
found  that  this  approach  produces  an  error  which  increases  progressively  as  the 
frequency  increases  at  a fixed  Mach  number.  At  higher  frequencies  it  appears 
that  the  ray  theory  approximation  is  likely  to  be  more  accurate  and  it  is  hoped 
to  publish  an  assessment  of  the  utility  of  this  method  later. 
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1 INTRODUCTION 

Experiments  have  shown  that  the  refraction  caused  by  vortices  shed  from 
lifting  surfaces  can  significantly  influence  the  noise  nuisance  of  an  aircraft* 
This  fact  has  led  to  vortex  refraction  being  the  subject  of  many  theoretical 
investigations,  in  order  to  understand  the  effect  better,  and  even  perhaps  to 
exploit  it. 

The  majority  of  published  analyses  are  based  on  ray  theory  whose  application 

2 

to  vortex  flow  has  been  welL  described  by  Broadbent  . Ray  theory  has  two  built-in 
restrictions  which  may  he  violated  in  practice.  First,  it  is  assumed  that  the 
acoustic  wavelength  is  small  compared  with  the  length  of  significant  variation  of 
the  underlying  flow  field.  Secondly  the  results  obtained  cannot  be  applied  in 
the  limit  of  the  extreme  far  field  if  the  predicted  acoustic  pressure  varies  by 
large  amounts  in  different  directions  due  to  such  effects  as  wing  shielding. 

Ray  theory  takes  no  account  of  the  smearing  out  of  discontinuities  like  shadow 
boundaries  which  takes  place  over  large  distances. 

The  first  restriction  is  not  serious  since  ray  theory  has  been  shown  to 
give  good  agreement  wilh  the  results  from  wind-tunnel  experiments  where  micro- 
phones normally  cannot  bo  placed  very  far  away  from  the  model.  However,  if  one 
wishes  to  predict  an  aircraft  noise  footprint,  for  example,  then  most  of  the 
regions  of  interest  on  the  ground  wiLl  probably  be  well  into  the  far  field. 

Rroadbont  also  mentioned  a second  theoretical  approach,  which  was  tried  as 
long  ago  as  1 59  ^ , based  on  the  assumption  that  the  maximum  Mach  number  M of 
the  vortex  flow  is  small  compared  to  one.  This  is  a reasonable  approximation  to 
make  for  an  aircraft  which  has  just  taken  off  or  is  about  to  land.  The  method 
consists  of  expressing  the  acoustic  pressure  as  an  asymptotic  power  series  in  the 
Mach  number.  In  practice  this  process  has  so  far  only  been  applied  to  first  order 
in  M , since  the  higher  order  terms  get  progressively  more  difficult  to 
calculate . 

Both  theoretical  approaches,  which  have  just  been  outlined,  pose  various 
questions  as  to  their  range  of  practical  application.  The  purpose  of  analysing 
in  detail  the  particular  example  described  in  Station  4 (Which  is  not  physically 
realistic,  but  is  comparatively  easy  to  solve)  was  to  assess  the  errors 
introduced  by  the1  different  approximations  that  can  be  made.  The  analysis  is 
greatly  simplified  by  the  use  of  a low  Mach  number  approximation  which  is 
uniformly  valid  for  all  frequencies.  This  contrasts  with  the  last  described 
method  of  approximation  which  clearly,  from  the  form  of  the  equations  to  bo 
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so  l veil  and  the  results  present  oil  in  iliis  Report,  is  not  uniformly  valiil  .is  t ho 

frequency  t i'tuls  to  intinitv.  The  results  trow  tin*  unitormly  valiil  approximat  ion 
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will  contain  an  error  ot  0(M~)  when  compared  with  the  exact  solution,  hut  the 
author  believes  that  this  error  is  always  insignificant  tor  the  range  ot  M 
considered  (from  0 1 1'  l'.  O.  Thus  the  results  described  in  section  t>,  for 
various  Mach  numbers  and  wave  numbers,  serve  as  a check  on  any  technique  used  to 
analyse  the  effect  ot  more  realistic  low  Mach  number  t low  fields  on  acmist  ie 
propagat ion. 

This  Report  only  contains  an  assessment  ot  the  asymptot  ic  power  series  in 
M for  acoustic  pressure,  but  it  is  hoped  to  produce  an  assessment  ot  ray  theory 
in  a future  pub  I teat  ion. 


THE  ACOUSTIC  EQUATIONS  IN  STEADY  TWO-DIMENSIONAL  HOMENTROPIC  VORTEX  FLOW 


In  what  follows,  the  adjective  ’basic’  refers  to  the  t low  in  the  absence 
of  any  acoustic  disturbance.  The  acoustic  equations  were  derived  in  Ket  •>  and  ate 
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where  t is  the  time,  K basic  densitv,  Q basic  velocity,  1’  basic  pressure, 
S basic  specific  entropy,  p is  the  acoust  ic  density,  q aceuat  ic  velocity, 
p acoustic  pressure,  and  s is  the  acoustic  specific  entropy.  It  viscous  and 
heat  conductivity  el  feet s are  ignored,  the  usual  equations  ot  fluid  mechanics 
show  that  the  basic  state  quantities  satisfy 
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Substituting  thoso  part  iv'itlai  torms  t't  t ho  basic  stato  ijuant  it  ton  into  equations 
l O to  (41  gives 
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« 0 , rhon  equal  ii'n  tin  shows  that  s is  constant  tor  oach  t luid 
tl  r 

particle  tor  all  time  and  thus  is  independent  ot  the  acoust  ic  dist  111  banco . 
implies  that  s must  bo  zero  tor  the  homont  topic  l low  considered.  Further 
equal  i i'ii s tSl  anil  HO. 
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Kqu.it  ions  (IS)  .mil  (I‘*1  may  now  In-  used  to  eliminate  v .inJ  w from  oqu.it  ions 
(IM  .nul  (171  anil  tlms  obtain  tin*  following  two  equations  wtiioh  couven iont  1 v 
.losoribe  acoustic  propagation  in  steady  homentropi.  vortex  tlow: 
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i [HI  U'OUST  1C  EQUATIONS  IN  FLUID  UNDKRGOIN(  SOUP  BODY  ROTATION 
I 1 In*  oxao  t equal  i ons 

Examinat ion  ot  oquat  ions  (.’01  and  (.'11  shows  that  they  have  their  simples 
form  when  V varies  I i near  I v with  r . tor  then  1)  and  the  operator  on  u i 
equal  ion  (’ll  heeome  independent  of  r . and  equations  (J03  and  y 1 3 can  bo 
easily  combined  to  give 


3 0 3 0 . . .... 

where  D is  now  . ♦ and  V , the  basic  velocity  disti ibut  ion 

3 t a 33  a 

takes  the  form  for  a solid  bodv  rotation.  Equation  (.'II  becomes 
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1 . 2 Separating  the  variables 

From  now  on  wo  will  be  dealing  with  an  acoustic  disturbance  which  is 
independent  of  z , and  has  a constant  frequency  w . Accordingly,  the  acoustic 
pressure  and  radial  velocity  component  can  bo  written  in  the  forms 
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Substituting  expression  (24)  into  equation  (22),  and  equating 


Fourier  components,  gives 
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where  C(a),  M and  R —■  . Similarly,  equation  (23)  gives 
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1.3  A uniformly  valid  approx i mat  ion  at  low  Mach  number 

Equation  (2b)  has  two  different  types  of  term  dependent  on  M . The  first 
is  contained  in  the  variation  of  C , which  can  be  quantified  using  the  equation 
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ignored  when  M ^ 1 , with  an  error  in  the  final  answer  which  is  truly  0(M*")  . 
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The  second  type  of  dependence  is  characterised  by  the  parameter  Mn  which  may  or 
may  nor  he  small,  depending  on  the  value  ot  n . It  is  therefore  important  to 
retain  all  the  terms  containing  Mn  if  serious  errors  are  to  be  avoided  even  at 
quite  small  Mach  numbers.  Kquation  (2b),  then,  will  be  approximated  by 
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which  is  a form  of  Bessel's  equation.  The  solution  of  this  equation  which  is 
finite  at  the  origin  is 
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where  A is  an  arbitrary  constant,  .1  is  the  Bessel  function  of  the  first 
n } n 


kind  and  nth  order,  and 


/ ■>  ’ / 

J(Ka  + Mn) “ - 2M“  1 + - 
n \ h 


— 

(a  + Mn  / 


The  advantage  ot  employing  the  above  approximation  is  clear  in  this  case, 

since  f can  then  be  evaluated  in  terms  of  a function  which  is  well  documented, 
n 

However,  the  author  believes  that  the  classi f icat ion  of  the  terms  containing  M 
into  two  types  is  important  in  general,  and  only  terms  which  are  uniformly 
small  may  safely  be  ignored  without  significant  error.  This  point  of  view  is 
reinforced  by  the  results  of  section  5 where  a cruder  approximation  process  is 
used  for  comparison. 

4 SCATTERING  BY  A REGION  OF  FLUID  IN  SOLID  BODY  ROTATION  OF  SOUND  GENERATED 
BY  A CON STANT  FREQUENCY  L INE  MONOPOLE 

4 . 1 The  mathemat ical  mode l 

The  particular  problem  which  we  are  going  to  investigate  in  detail,  using 
the  approximation  developed  in  the  last  section,  is  illustrated  in  Fig  I.  It 
consists  of  a region  of  fluid  in  solid  body  rotation  bounded  by  the  cylinder 
r,  = 1,  with  an  infinite  region  of  stagnant  fluid  outside  containing  a monopole  of 
acoustic  pressure  on  the  line  given  by  = 5 and  d = n . The  possibility  of 
instabilities  developing  on  the  vortex  sheet  separating  the  two  regions  will  be 
ignored  in  this  analysis. 

The  approximation  described  in  section  1.3  is  uniformly  valid  since  the 
region  of  solid  body  rotation  is  bounded.  Thus  the  acoustic  pressure  is  given  by 
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when  £ 1.  In  the  stagnant  fluid,  the  acoustic  pressure  is  the  sum  of  two  out 

going  waves,  one  due  to  the  monopole,  the  other  to  scattering  by  the  rotating 
fluid.  The  first  contribution  is  the  incident  wave  which  represents  the  total 
acoustic  pressure  in  the  limit  M ► 0 . This  is  given  by 


V>1'X>^(Ka/^  * C0  + “U0  cos  fl)  * 


where  V.  is  a eiven  const  ant  and  H is  the  liankol  timet  ion  oi  the  second 
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kind  and  nth  order.  Kquat ion  (33)  can  be  rewritten  as 


Pi  - 1 


I 


f-  l)nH(J)(Kai.  ,)J  (KaOeind 
n 0 n 


(i<i.0>.  (34) 


by  using  Graf's  addition  theorem.  (See  Ref  1,  result  ll.|,7s.)  Then  the  total 
aeoustic  pressure,  in  the  region  I » J » Cq  , is  given  by  an  infinite  series  of 


the  form 
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where  B is  an  arbitrary  constant, 
n 

Since  the  expansion  in  equation  (34)  is  required  for  the  application  of 
boundary  conditions  at  ■ 1 so  that.  A and  B can  be  determined,  we  have 
chosen  the  form  which  is  valid  in  the  region  e,  < f.„  . This  restriction  does 
not  affect  the  subsequent  evaluation  of  the  far-field  acoustic  pressure  described 
in  section  4.3. 

4.3  t'ond i t ions  on  the  vortex  sheet  hound i ng  t he  solid  body  rotation 

There  are  two  boundary  conditions  which  are  normally  applied  across  a 
vortex  sheet.  They  are: 

(i)  continuity  of  pressure; 

(ii)  continuity  ot  the  velocity  component  normal  to  the  sheet. 
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The  first  condition  ensures  that  the  acceleration  of  fluid  particles  in  the 
sheet  is  finite,  and  the  second  is  a consequence  of  the  principle  of  conserva- 
tion ot  mass.  Assuming  that  the  basic  pressure  has  been  made  continuous, 
condition  (i)  implies  that  p is  continuous  across  the  sheet  and,  in  particular, 
the  Fourier  components  in  equations  (Id)  and  (35)  are  equal  at  ( * 1,  :< 
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Condition  (ii)  is  more  complicated  to  apply.  Suppose  the  surface  of  the 
vortex  sheet  is  described  by  = 0.  The  form  of  I should  show  that  the 
aeoust ic  disturbance  linearly  perturbs  t He  sheet  away  from  the  circle  r = a 
which  the  sheet  would  occupy  if  the  lino  source  were  absent.  Thus  I may  be 
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For  a point  r on  the  ; aet  t , 
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Condition  (ii)  implies  that 
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where  £>-1  + means  that  the  unit  circle  is  approached  from  outside,  and 
f.*d-  means  that  the  unit  circle  is  approached  from  inside.  After  rejecting  the 
products  of  linear  perturbations  and  using  the  definition  of  I , these  equations 
simplify  to 


L i ra  u 

C-I  + 


> 

dr 

dt 


vi: 


Lim  u 
4*1- 


iV 


0 


a 


iwt 


I 


n*-*» 


. inO 
no  e 
n 


1 1 


(41) 


Thus,  combining  equations  (39)  and  (41),  using  equation  (25),  and  equating 
Fourier  components,  we  obtain 
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which  leads  to 
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Now,  equation  (37)  applies  outside  the  solid  body  rotation  if  M is  made 


.•cro,  and  so  it  can  be  used  to  expn 
equation  (43)  to  give 
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(44) 


Finally,  substituting  the  forms  of  f i nip  1 i e«J  by  equations  (32)  and  (35),  we 

obtain 

R ( a ) A 


Ka  (Ka  + Mn)ci  J ' (a  ) + 2MnJ  (ci  )> 
n 1 n n n n n j 

HU>  (Kaf.n)(Ka  + Mn)  |{Ka  + Mn)2  - 4M~MJ'(Ka)  + B H(2)'(Ka)> 
n 0 I j I n n n J 


= P. (-  1) 
i 


F.quations  (3b)  and  (45)  together  determine  the  values  of  A and  B 

n n 


(45) 


4 . 3 Calculating  the  effect  on  the  far-field  acoustic  pressure 

We  are  now  in  a position  to  calculate  the  acoustic  pressure.  In  order  to 
keep  the  number  of  calculations  within  reasonable  bounds,  numerical  results  were 
obtained  only  for  the  far  field  where  the  azimuthal  variation  of  p is 
independent  of  £ . This  can  be  seen  by  taking  the  limit  as  £ -*■  » of  the 
right-hand  side  of  equation  (33)  and  using  the  principal  asymptotic  form  of 
(described  by  result  9.2.4  of  Ref  3)  so  as  to  obtain 
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2_  i {a)t-Ka  (f.+f.Q  cos  0)  + (tt/4)  } 

P " i'vj  uKaf,  e 


I,  iKaf.n  cos  0 
x <1  + e ») 


(-l)nBnH<2)(Kat0)ei"9  * 0 jdf  • J 


1 'o  Ka'-o 


The  effect  of  the  solid  body  rotation  on  the  far-field  azimuthal  variation 
of  p was  assessed  by  evaluating  the  ratio  of  the  amplitudes  of  p and  , 

Cc  the  far-field  directivity  defined  by 


10  I op,  -2-  - io  log  1 + eiKa^COS0 

1 0 p . • 10 

1 


(-  i)n  B n(2)  (KaCn)elnL* 
n n 0 


Eliminating  A from  equations  (.36)  and  (45)  gives 


3 J (Ka)  - Y J’ (Ka) 
g n_  n n_  _n_  _ _ 

Y H(2)  ’ (Ka)  - B (Ka) 

n n n n 


where 


B - Ka|(Ka  + Mn)a  J'  (a  ) + 2MnJ 
n l n n n 


Y = (Ka  + Mn)  l(Ka  + Mn)'  - AM2 


id  (a  )} 
n n I 

!}d  (a  ) 
n n 


The  infinite  series  in  equation  (47)  has  similar  convergence  properties  to  that 

which  would  be  obtained  if  the  rotating  region  were  replaced  by  a hard  circular 

cylinder.  For,  as  n -►  »«>,  y /p.  ► ■*  and 

n n 


(-  i)  nB  H(2) (Kaf  )eln6 
n n 0 


J ' (Ka)  ...  . . 

/ • \n  n „(2 ) , mG 

- (-  O -Tvv H (KaCn)e 

II  2 ' (Ka)  " ° 

n 


which  is  the  exact  form  of  the  nth  term  in  the  series  representing  the  far- field 
acoustic  pressure  scattered  when  -^2  = o on  r = a . Since  the  latter  series  is 
known  to  converge,  it  follows  that  the  corresponding  series  in  equation  (47)  also 
converges.  There  is  a further  similarity  between  the  two  series  in  that,  to 
obtain  a required  accuracy,  progressively  more  and  more  terms  have  to  be  summed 
as  Ka  increases. 


/c  J. 
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1 1 


I'lit*  series  in  ei|uat  ion  l-WI  was  approx  imat  ed  l>v  C ho  pail  ial  -.urn  I i om 
n - m to  rm  , wlie if  m l-  1)1  was  olios, m so  that  the  terms  omitteil  were 
smaller  in  magnitude  than  the  lomnlinjA  error  ol  the  eomputei  upon  whieh  the 
ea  leu  tat  ions  were  pet  termed.  litis  |>roeetlnre  was  just  it  ie,l  hv  ineroasiuy.  m ami 
showing  that  the  ettert  was  indeed  nee.  1 i ,•  t h 1 e 1 1>  the  order  ot  aeeuraev  ot  the 
eomputer  in  single  preeision  ar  i t hmet  i e . the  results  are  dismissed  in  seetion  t< . 


I'ltl''  KVAl.UAT  ION  OF  Fill  U’Ol'STIC  I’KFKSt'KI  AS  AN  ASYMPTOTIC  I'OWl'K  SIKHS  IN 
MACII  Nil  Mill  K 


. I A^ip  1 ieat  i on  ol  the  method  to  vortex  ! I ow  in  peueial 


We  now  return  to  eijuat  ions  l-’O)  and  l I ) , and  st>lv«‘  them,  when  M I 
hv  the  see, nut  theoretieal  appioaeh  outlined  in  the  i ut  rodue  t ion . I.et 


Fill  V'V,  where  V, 
0 0 


max  V , and  lei  C,  he  a e ha  i ae  t <•  l i st  i c value  ot  the 
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speed  ot  sound.  llie  detiuit  ions  ot  \'  , and  C , aie  t ompnt  it»le  with  those 

0 0' 


adopteil  in  seetion  ' provided  F - toi  t h,'  » l i , l hodv  rotation.  hsinp.  these 
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where  A ’*  > M 

V1  , .10 


t',|uat  ions  i'M  and  t,li)  imply 


,1 


and 


'U'/c 
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K / .V\  V 


c V ■'  i,  j 


vd 


S i"  i 


so  that  K and  C are,  tormal It  at  least,  eonstant  to  OiMl 


In  this  met  hod,  it  is  assumed  that  p and  u ran  he  expressed  as 
asvmptot  it'  power  set  ies  in  M , 
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wlH'ii'  p()  ami  u(  an>  independent  of  M . Stibaf  i tut  i nK  t hose  series  into 
equations  (•>.’)  an.l  (SI),  not  in*  ,'i|uat  ions  (54)  ami  (55),  an, I equation  powers  ot 
M , w,<  timl  that  til,'  first  two  tonus  are  ^ i von  l>v 
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Mto  i',|nat  ion  lor  p^  , whiolt  is  t ho  inoiilont  aoonst  io  pressure,  is  t ho 
or, 1 1 nary  wavo  oquat  ion  p.ovorn  i n«  aoonst  io  propagation  in  a stationary  homoo.onoous 
mo, limn.  :V  t ho  hasio  stato  ol  t ho  tlnid  when  M - 0.  Tho  first  order  olloot  ot 
t ho  Mow,  in  this  mo  do  l , is  to  intvoduoo  a sontoo  d i st  r ihnt  ion  (tho  rip.lit  hand 


I0t' 


oni*/  h 


side  ot  equation  IbL'))  w'aiiti  is  tin  result  ol  interaction  between  t he  given 
incident  acoust  ic  pressure  and  tli.it  I low.  The  effort  ol  the  vortex  can  bo 
expressed  in  a similar  fashion  to  equal  ion  (47),  that  is 


10  log  „ 
' 10 


l’| 

It)  M lop.,  e " Re  + 0(M* ) 

10  P0 


where  Ke  denotes  the  real  part. 

5*2  The  method  appl ied  ti>  the  problem  ol  section  4 

The  bounded  region  of  solid  body  rotation  in  the  analysis  ot  section  ■> 
implies  a non-dimensional  azimuthal  velocity  component  given  by 


K - 1 11  (a  - r ) 

a 


(04) 


where  the  Heaviside  function  H(x)  is  defined  by  ll(x)  m l for  x **  0 and 
H(x)  « t'  for  x v 0 ; here  it  arises  as  a consequence  of  the  vortex  sheet  at 
r " a . The  incident  acoustic  pressure  is  given  by  equal  ion  (3  1),  and  so 
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0 


T.e 
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i -c  t 


(bS) 


which  is  a solution  of  equation  (hi).  Substituting  equation  (b4)  into  equation 
(b2),  we  obtain 


1 JL 

r 3r 
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'S  > 

r “ 3 0 ‘ 
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Ka 


K'  Ufa 


r) 


+ .Hr 


(ho) 


where  use  has  been  made  of  the  fact  that  Pj  is  independent  of  z and  has 
frequency  w . The  factor  Mr  - a)  , where  5 is  Dirac's  delta  function,  arises 
from  differentiating  H(a  - r)  . 


The  solution  of  equation  (hh)  can  be  found  by  convolving  the  source  term 
with  the  two-d  imetis  i ona  l free  field  (Iroeu's  function  for  Helmholtz's  equation 
(appropriate  to  outgoing  waves'),  namely 
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r . =0  0 . =-n 
1 l 


K"H(a  - r . ) + 6 (r 


. - a)  (~ — ] ■h[  ^(k/t^  + r2  + 2r.r  cos  0.)dr.d6. 

l ^3r.  a J 0 l 0 l 0 \)  l l 


and  Lf  use  is  made  of  Graf's  addition  theorem  as  in  equation  (34),  then 
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K2  / r.J2(Kr.)dr.  + J (Ka)  (KaJ'(Ka)  - J (Ka) 
J l n i i n [ n n 


The  integration  with  respect  to  r.  can  be  performed  using  results  I 1.3. 31, 
9.1.5  and  9.1.27  of  Ref  3 to  give 
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Finally,  using  equation  (63)  and  the  principal  asymptotic  form  of  ' , we  can 

show  that  in  the  far  field 
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where  Im  denotes  the  imaginary  part. 

There  are  a number  of  interchanges  of  limiting  processes  in  going  from 
equation  (67)  to  (72).  They  can  all  be  shown  to  be  justified  since  the  region  of 
integration  is  effectively  finite  and  the  series  are  all  uniformly  convergent. 

It  can  also  be  shown  that  the  form  of  the  nth  term  in  equation  (72)  is  the  same 
as  that  which  would  be  obtained  from  equation  (47)  by  writing  it  correct  to 
0(M)  . This  provides  a check  on  the  working  used  to  derive  both  equations. 

The  most  interesting  feature  of  the  final  expression  is  that  it  is  an  odd 
function  of  0 . This  is  a general  property  of  10  log.J-^-1  , correct  to  0(M)  , 


is  an  even  function  of  0 


This  can  be  deduced  from  equation  (62), 

3Po 


since  its  source  term,  which  is  a function  of  -rr-  , must  be  odd,  and  its  operator 

. . . 36  2 
on  P|  preserves  parity.  This  feature  has  also  been  demonstrated  by  Broadbent  , 

starting  from  the  expressions  obtained  in  the  ray  theory  limit. 
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t<  DESl  RW'l  lON  AND  DISCUSSION  OK  KKSUI  TS 

t’ . I The  pros  out  at  ion  ot  results 

the  tiLtin  body  ot  results  presented  in  Figs  2 to  !t>  .ire  calculated  from 
equation  (47).  The  values  of  tin-  dependent  variable  are  shown  only  for  the  range 
ot  0 from  -4S'  to  ♦45*’  where  the  vortex  has  the  most  significant  effect  on  the 
acoustic  pressure.  In  this  particular  example  there  is  also  a significant  effect 
in  the  region  on  the  opposite  side  of  the  solid  body  rotation  because  t !i»-  vortex 
sheet  at  the  edge  ot  the  rotation  causes  back  scattering  in  a manner  analogous  to 
a solid  circul.it  cylinder.  However,  this  region  is  ignored  since  the  hack 
scattering  would  not  be  present  in  more  practical  cases  where  the  basic  velocity 
t ie Id  i . coni inuous . 


A further  restriction  on  the  results  presented  here  is  that  they  are  all 
for  a single  value  of  (|  (f,()  1 2)  so  as  to  study  in  detail  the  variations 

brought  about  by  changes  in  Mach  number  and  wave  number.  Some  cal eul at  ions , 
which  were  done  for  higher  values  of  t.()  , indicate  that  the  directivity  pattern, 
tor  fixed  values  of  M and  Ko  , is  reduced  both  in  amplitude  and  angular  spread 
as  the  distance  of  the  source  from  the  centre  of  rotation  increases. 


The  graphs  of  far  field  directivity  arc  grouped  in  three  different  ways 
depending  on  the  range  ot  Mach  number  and  wave  number.  Figs  l?  to  7 art'  designed 
to  show  how  the  directivity  changes  as  M -*•  0 at  fixed  Ka  . Since  equation  (72) 


is  independent  of  M , the  two  curves  calculated  from 


shows  that  him  logl-^— 

M81  M |P0 

equation  (47),  at  each  fixed  wave  number,  are  scaled  with  Mach  number  so  that 
they  should  he  close  to  the  third  curve  which  is  calculat  ed  from  equal  ion  (72) 
for  a Mach  number  of  0.1.  The  latter  value  of  M is  chosen  for  convenience, 
since  then  ' log.  . -E- 

M 10  P() 


equals  the  far-field  directivity  defined  by  equation  (47). 


Figs  8 to  10  and  II  to  lb  both  have  the  same  dependent  variable  calculated 
t rom  equal  ion  (47),  but  the  former  group  emphasises  dependence  of  far-field 
directivity  on  wave  number  whereas  the  latter  group  compares  directivities  at 
fixed  Ka  , again,  for  the  higher  Mach  numbers.  At  these  values  of  M , the  low 
Mach  number  approximat ion  made  in  sect  ion  1 is  probably  beginning  to  break  down 


6 . 2 1'he  variation  of  far-tield  directivity  with  Mach  number  and  wave  number 


I 


I 

i 


The  dependence  of  far-field  directivity  on  Mach  number  and  wave  number  is 
shown  to  he  quite  complicated.  However,  there  are  some  trends  which  can  he  g 

identified  from  Figs  2 to  16.  As  M » II  with  Ka  fixed,  the  directivity  tends  S? 

a 

to  the  odd  function  of  0 predicted  hv  the  analysis  ot  section  8 correct  t o 0(M)  100  K 


I 'I 


but,  as  M increases  away  from  zero,  the  directivity  becomes  less  l i Wo  an  odd 
(unction.  For  the  lower  values  ot  Ka  , lor  wliieb  the  directivity  pattern 
remains  relatively  simple,  this  divergence  may  be  charact or a sed  by  the  behaviour 


ot  the  three  parameters  0 I 

1 0 max 


t lit'  smallest  positive  value  ot  0 where  10  1 oy 


and  1 which  are  defined  respectively  bv 

min  , 

L\ 


10 


value  of  It'  log 


10 


P; 


is  /ert>,  t lie  hum  mum 


, and  tin'  minimum  value  ot  10  log 


10 


riit' se  pai  amt* tors 


are  plotted  against  Mach  number,  tor  ka  - S , |0  and  1 k , in  Fig  17.  I'll  is  shows 


that  as  M increases. 


0 


increases  away  t torn  the  value  zero  predicted  by  the 


0(M)  theory  ot  section  S,  and  that  l 


max 


is  increasingly  overpred  i c t etl  and 


I I increas i ugly  under -predict ed  bv  the  Of  Ml  theory. 

min'  c 


For  all  the  values  ot  M and  ka  covered  by  Fig  17,  the  directivity  can 
be  simply  described  as  giving  a region  ot  decreased  aconst ic  pressure  t oi  0 «•  0 ^ 
fol lowed  by  a region  ot  increased  aconst  ic  pressure  tor  0 " 0.  . The  value  ot 
M , upto  which  this  simple  picture  holds,  decreases  as  ka  increases  and  the 


direct i v i t 

i y c 

'voices  into 

an  intertereuce  pattern  wtii 

eh  is 

probably  bet t ei 

ill' Si'  1 i bed 

by 

ray  theory. 

l’he  numbers  ot  maxima  and 

mini ma 

, in  the  region 

- to"  < o' 

> 

< 

10" . are  pi, 

,'t  t ed  against  wave  number  at 

t i xed 

Mach  number  in 

V i\\ 

I'll  is  clcarlv  shows  that  as  ka  increases,  the  numbers  increase  and  thus  t he 
angular  distance  between  direct  ions  ot  constructive  and  destructive  intertereuce 
decreases,  an  important  property  ot  int ert et once.  It  is  hoped  to  give  a detailed 
explanation  ot  how  this  intertereuce  pattern  arises,  as  well  as  an  assessment  ot 
the  accuracy  of  rav  theory  in  this  example,  in  a tuture  publicat ion. 


The  overall  implication  ot  these  trends  is  that  the  asymptotic  expansion 
in  powers  ot  M is  not  uniformly  valid  as  ka  ' . The  reason  toi  this  can  be 

readily  seen  from  the  remarks  made  on  the  convergence  ot  the  infinite  series  in 
equation  (.  •»  7 > . As  ka  increases,  more  terms  have  to  be  included  to  achieve  a 
given  accuracy.  This  implies  that  the  average  value  ot  |Mn|,  amongst  the 
significant  terms,  increases  and  it  reaches  a point  where  the  effect  of  this 
parameter  can  no  longer  be  regarded  as  small  enough  for  the  asymptotic  power 
series  approach  to  be  adequate. 


More  generally,  ray  theory  shows  that  it  is 
the  phase  and  amplitude  ot  the  acoustic  pressure, 
the  method  of  section  '■>  is  t o adopt  an  asymptotic 
both  amplitude  and  phase  when  ka  ^ I . 


necessary  to  distinguish  between 
Vims  an  obvious  improvement  to 
expansion  m powers  ot  M tor 


10(1 


.'0 


t> . ' Assessment  .'1  ! lie  utility  ot  the  is  vmpt  ot  i c power  si'rirs  in  M lor 

■nous  t i e pres.su  ft* 

In  n'ot  ion  S,  t hi1  .is  vmpt  ut  i i expansion  was  taken  only  to  t)(M)  . In  the 

example  used  in  the  present  Repot t however,  that  ordei  of  accuracy  was  achieved 

with  about  the  sans-  amount  ot  computational  ottort  as  was  required  for  the 

inherently  mote  accutate  method  ot  analysis  described  in  sections  1 and  4.  To 

2 

take  the  asymptotic  expansion  a stage  further  to  would  require  at  least 

twice  the  ottort  needed  tor  the  other  approach.  Although  this  may  not  be  a 
general Iv  applicable  criticism,  in  that  the  t unctions  derived  in  section  3 are 
particularly  simple  tor  the  solid  body  total  ion,  the  results  have  shown  that  the 
asymptotic  expansion  is  not  uniformly  valid  as  ka  • •»  , and  so  this  mo  t hod  will 
not  give  necessarily  useful  estimates  ot  the  vortex  retraction  eltect  in 
prac t icu l cases . 

It  seems  to  the  author  that  the  best  approj.  h to  the  general  problem  of 
predicting  the  effect  of  a vortex,  a l indeed  the  pt  opag.it  ion  ettects  ot  all  low 
Mach  number  flows,  is  to  retain  the  0(M>  terms  in  the  operator  on  acoustic 
pressure  and  then  attempt  to  solve  the  resulting  approximate  acoustic  equation  as 
accurately  as  possible.  An  example  ot  the  sucecstul  tpp  i cat  ’ on  ot  this 
approach,  in  the  case  ot  potent  i at  t I.  , is  Jcsct  i cd  in  lii  I s. 

7 t own  s ions 

Hie  scattering  ot  sound  by  a rigion  ot  tluid  in  solid  bodv  rotation  at  low 
Mach  number  is  an  easily  analysed  i x.cnple  which  . i.i  he  used  as  i test  case  fot 
assessing  the  accuracy  ot  various  ap;  oxi"iat  e met  hods  tot  c.  ilculat  ing.  the  effect 
on  aeon  st  ic  propagat  i on  ot  mor  • pt  illicit  vort  \ ’ 1 ou  ■ Idn  main  body  ot  results 
presented  in  this  Report  are  obi  lined  bv  solvin’,  the  acoustic  cqu.it  ion  after 
ignoring  certain  terms  which  arc  OfM'  1 and  arc  believed  to  h.iv«  a uniformly 
small  cited  at  least  up  to  the  max i mum  Mach  number  considered  ot  0.  1.  for  this 
reason,  the  method  is  considered  to  be  suitable  fot  use  as  a check  against  which 
other  approxim.it  e methods  can  he  tested. 

One  such  method  is  that  of  an  asymptotic  expansion  m powets  ot  M . It 
has  been  shown  that  such  an  expansion  is  not  unitormlv  valid  lot  large  ka  and, 
when  tested  against  the  uniformly  valid  results,  it  i s t ouml  to  be  seriously  in 
error  for  ka  * IS  even  at  a Mach  number  ot  0.1. 

for  large  wave  numbers  the  directivity  ot  acoustic  pressure  takes  on  the 
appearance  ot  an  interference  pattern.  This  implies  that  rav  theory  will  provide 
a better  descript  ion  at  these  wave  numbers  and  the  next  task  will  be  to  assess  ^ 

the  accuracy  of  this  approximation. 
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